Tuning the Resonance in High Temperature Superconducting Terahertz Metamaterials 
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In this Letter we present resonance properties in terahertz metamaterials consisting of a split-ring resonator 
array made from high temperature superconducting films. By varying the temperature, we observed efficient 
metamaterial resonance switching and frequency tuning with some features not revealed before. The results 
were well reproduced by numerical simulations of metamaterial resonance using the experimentally measured 
complex conductivity of the superconducting film. We developed a theoretical model that explains the tun- 
ing features, which takes into account the resistive resonance damping and additional split-ring inductance 
contributed from both the real and imaginary parts of the temperature-dependent complex conductivity. The 
theoretical model further predicted more efficient resonance switching and frequency shifting in metamaterials 
consisting of a thinner superconducting split-ring resonator array, which were also verified in experiments. 
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Metamaterials consisting of metallic elements have enabled 
a structurally scalable electrical and/or magnetic resonant re- 
sponse, from which exotic electromagnetic phenomena absent 
in natural materials have been observed [1]. Metals provide 
high conductivity that is necessary to realize strong electri- 
cal/magnetic metamaterial response EJO. Metals, however, 
play a negligible role in active/dynamical metamaterial reso- 
nance switching and/or frequency tuning, which has been typ- 
ically accomplished through the integration of metamaterials 
with other natural materials (e.g. semiconductors) or devices, 
and by the application of external stimuli |4" MT2| . It is essen- 
tially the modification of the metamaterial embedded environ- 
ment that contributes to such previously observed functional- 
ities. 

Recently, there has been increasing interest in supercon- 
ducting metamaterials towards loss reduction Ifl3l - l20l . Signif- 
icant Joule losses have often prevented resonant metal meta- 
materials from achieving proposed applications, particularly 
in the optical frequency range. At low temperatures, su- 
perconducting materials possess superior conductivity than 
metals at frequencies up to terahertz (THz), and therefore 
it is expected that superconducting metamaterials will have 
a lower loss than metal metamaterials. More interestingly, 
superconductors exhibit tunable complex conductivity over a 
wide range of values, through variation of temperature and ap- 
plication of photoexcitation, electrical currents and magnetic 
fields. Therefore, we would expect correspondingly tunable 
metamaterials, which originate from the superconducting ma- 
terials composing the metamaterial, in contrast to tuning the 
metamaterial environment. 

In fact, superconducting metamaterials have enabled dia- 
magnetic response at very low frequencies, which may en- 
able screening of static magnetic fields ET1 |22l . In the mi- 
crowave frequency range (^10 GHz), left-handed supercon- 
ducting transmission lines have been introduced and their tun- 
ability has been realized with electrical currents or temper- 
ature iTPJl [141 . Negative index metamaterials comprised of 
niobium wires and split-ring resonators exhibited red-shifting 



in resonance frequency (~ 0.6%) when the temperature was 
increased approaching the transition temperature T c ifTTl . Fur- 
ther experimental work has demonstrated tunability through 
application of external dc or rf magnetic fields ifTBI |20l . 
At THz frequencies, low temperature superconductors may 
be not suitable for metamaterial applications, since with a 
smaller superconducting gap, Cooper pairs may be excited 
and broken by THz photons, and therefore high temperature 
superconductors (HTS) should be employed. In this Letter, 
we present THz metamaterials based on electric split-ring 
resonators (SRRs) made from epitaxial YBa2Cu3C<7_8 HTS 
films. These metameterials exhibit temperature-dependent 
resonance strength and frequency, which reveal some inter- 
esting tuning features not previously observed. Finite-element 
numerical simulations and theoretical modeling are performed 
to understand the underlying tuning mechanism. 

The epitaxial YBa 2 Cu 3 07_s (YBCO) films with 8 = 0.05 
were prepared using pulsed laser deposition on 500 /jm thick 
(100) LaAlC<3 (LAO) substrates. The transition temperature 
was measured to be T c = 90 K. Square arrays of electric SRRs, 
with the unit cell shown in the inset to Fig. [T|a), were fabri- 
cated using conventional photolithographic methods and wet 
chemical etching of the YBCO films. The YBCO SRRs have 
a thickness of d = 180 nm or 50 nm, outer dimensions of 
/ = 36 jum, line width of w — 4 /jm and gap size of g = 4 fjm, 
and the arrays have the periodicity of p = 46 pm. Terahertz 
time-domain spectroscopy (THz-TDS) incorporated with a 
continuous flow liquid helium cryostat was used to charac- 
terize the YBCO films and metamaterials. Under normal inci- 
dence, the THz transmission spectra were measured as a func- 
tion of temperature, using an LAO substrate as the reference. 

We focus our attention on the metamaterial fundamen- 
tal resonance (the so-called LC resonance) resulting from 
the circulating currents excited by the incident THz electric 
field ll23l . In Fig.[TJa) we show the THz transmission ampli- 
tude spectra for the 180 nm thick YBCO metamaterial sam- 
ple at various temperatures. At temperatures far below T c , 
e.g. 10 K, the metamaterial exhibits the strongest resonance, 
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FIG. 1: (color online), (a) THz Transmission amplitude spectra of 
the 180 nm thick YBCO metamaterial at various temperatures, (b) 
Transmission minimum and (c) corresponding resonance frequency 
as functions of temperature, from experiments, numerical simula- 
tions, and theoretical calculations. Inset to (a) illustrates a micro- 
scopic image of a single YBCO SRR, where the lighter colored area 
is YBCO. Inset to (b) shows the real and imaginary parts of the com- 
plex conductivity at 0.6 THz of an unpatterned 180 nm thick YBCO 
film. 



as indicated by the sharp THz transmission dip with a min- 
imal transmission amplitude of 0.045 at 0.613 THz. This 
strong resonance is almost the same as in a metamaterial sam- 
ple where the YBCO SRRs was replaced by gold SRRs with 
the same thickness and at the same temperature. As the tem- 
perature increases, the resonance strength decreases, as seen 
by the broadening and reduction in amplitude of the transmis- 
sion dip. The resonance frequency experiences a red-shifting, 
which reaches the lowest value of 0.564 THz near 80 K, re- 
sulting in a frequency tuning of 8%. As the temperature fur- 



ther increases, the resonance strength continues to decrease, 
but the resonance frequency, on the other hand, shifts back 
to higher frequencies. The temperature-dependent transmis- 
sion minimum and the corresponding resonance frequency are 
plotted in Figs.[TJb) and[TJc), respectively. The results show 
that, at temperatures near 80 K, the transition of resonance 
strength is fastest and the resonance frequency exhibits a dip, 
not observed in previous work 1 17-20 1. We can exclude the 
LAO substrate as contributing to the metamaterial resonance 
tuning, because the features in the temperature-dependent res- 
onance in the YBCO metamaterial (see Fig. [TJ were not ob- 
served in a metamaterial sample where the YBCO was re- 
placed by gold. In that gold SRR metamaterial sample, the 
resonance frequency shifting was imperceptible, and the res- 
onance strength only slightly decreased with increasing tem- 
peratures. Additionally, through THz-TDS measurements, it 
turns out that the dielectric constant of the LAO substrate only 
exhibits a weak dependence on the temperature. Therefore, it 
is the temperature-dependent properties of YBCO film that are 
responsible for the observed metamaterial resonance tuning. 

The complex conductivity of YBCO film can be expressed 
using the well-known two-fluid model |24|: 
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where /„ and / s are fractions of normal (quasiparticle) and 
superconducting (superfluid) carriers, respectively, with f n + 
f s = 1, n is the carrier density, m* is the carrier effective mass, 
and x is the quasiparticle relaxation time. The real and imagi- 
nary parts of the complex conductivity are then: 
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Using THz-TDS we experimentally measured the conduc- 
tivity of an unpatterned 180 nm thick YBCO plain film. The 
resultant real and imaginary parts of the complex conductivity 
at 0.6 THz are plotted as functions of temperature in the inset 
to Fig. |TJb). The real conductivity [Eq. Q], which derives 
from the Drude response of quasiparticles in Eq. ([T]), slowly 
increases when temperature decreases across T c to about 70 K. 
It starts to decrease below 70 K, but not significantly, over 
the temperatures we measured down to 50 K. In this temper- 
ature range, cox <C 1 at the resonance frequency (~ 0.6 THz), 
the decreasing f n (T) may be compensated by the increasing 
quasiparticle scattering time x |25|. At temperatures above T c , 
the imaginary conductivity is derived from the Drude response 
and is very small, since f s (T > T c ) = 0. As the temperature 
decreases below T c , the second term in Eq. <[3j from the su- 
perfluid Cooper pair state becomes non-zero and results in the 
rapidly increasing imaginary conductivity, exceeding the real 
conductivity below 80 K. Using these experimental values of 
the YBCO complex conductivity at 0.6 THz, the metamaterial 
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resonant response was simulated using commercially avail- 
able finite-element simulation codes from COMSOL Multi- 
physics. The simulated transmission minimum and the corre- 
sponding frequency are plotted as functions of temperature in 
Figs.[TJb) and[TJc), respectively, reproducing the experimental 
results. 

The measured real conductivity of the YBCO film reveals 
less than 20% change over the temperature range from 60 K to 
90 K, where the resonance strength experiences a fast change 
and the resonance shifts its frequency. This variation of real 
conductivity cannot solely cause the observed large metama- 
terial resonance switching and frequency tuning. Both the real 
and imaginary parts of the complex conductivity have to be 
considered for the metamaterial resonance. The imaginary 
conductivity, which is due to the superfluid carriers and causes 
no loss, becomes dominant at low temperatures, and it is re- 
sponsible for the enhancement in resonance strength. 

The resonance frequency is determined by the effective 
capacitance, C, inductance, L, and resistance, R, in the 
SRRs [1251 : 



can be further simplified as: 
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It has been shown that the kinetic inductance, which repre- 
sents the kinetic energy storage in free electrons in metals, or 
Cooper pairs in superconductors, plays an important role in 
determining the metamaterial resonance frequency ifTTl 1271 . 
This effect underpins the red-shifting of the resonance fre- 
quency in niobium metal superconducting metamaterials op- 
erating near 10 GHz as the temperature increases and ap- 
proaches T c |[m . However, the back blue-shifting, shown in 
Fig. [TJc) between ^80 K and T c , was not observed in nio- 
bium, and the model proposed in that work would not explain 
this effect when only the superfluid state (i.e. imaginary con- 
ductivity) was considered 1 1 7|. 

Here we consider a more general situation that the SRRs 
are fabricated from a conducting film (YBCO film in our case) 
with a complex conductivity 6 and thickness d. Such an un- 
patterned plain film can be modeled as a lumped impedance 
in an equivalent transmission line. By equating the (multi- 
ple) reflections or transmissions from the film and the trans- 
mission line model, this complex surface impedance (or sheet 
impedance with units of il/square) of the unpatterned film can 
be derived as: 



Zs = Rs — iXs = Zq 
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where the tildes over the variables indicate complex values, 
Zo = 377 £2 is the vacuum intrinsic impedance, «3 = -^/Elao = 
4.8 is the LAO substrate refractive index, hi = J /6/eoCO is the 
complex refractive index of the film, and p = «2Co/co is the 
complex propagation constant where co is the light velocity 
in vacuum. Both «2 and p can be calculated from the experi- 
mental complex conductivity near the metamaterial resonance 
frequency. When ns -C \ni\, which is valid in our case, Eq. <j5j 
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From Eqs. Q and d§J, it is obvious that both the finite real 
and imaginary parts of the film refractive index H2, and there- 
fore the finite real and imaginary parts of the complex con- 
ductivity 6, contribute to the film surface resistance Rs and 
reactance Xg,. They are plotted in Fig.|2jfor the 180 nm thick 
YBCO film as functions of temperature, and are also calcu- 
lated for a 50 nm thick YBCO film assuming its complex con- 
ductivity does not depend on the film thickness. 

The YBCO SRR array resistance R (SRR reactance X is 
zero at resonance) can be obtained by considering the nonuni- 
form distribution of currents in a unit cell |28"1 : R = [(A — 
g)/w]Rs, where A — 64 /jm is the median circumference of 
the (small) current loop. When temperature increases, the in- 
creasing SRR resistance R accounts for the resonance damp- 
ing and therefore the increasing transmission minimum. If we 
model the SRR array as a lumped resistor R in the transmis- 
sion line 1291 , we can calculate the resonance transmission as 
a function of temperature, which is plotted in Fig.[TJb) and sat- 
isfyingly reproduces the experimental and simulated results. 

In order to correctly interpret the temperature-dependent 
resonance frequency shifting, additional inductance in SRRs 
has to be taken into account besides the geometric inductance 
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FIG. 2: (color online). Complex surface impedance of the 180 nm 
thick unpatterned YBCO superconducting film calculated using the 
experimental complex conductivity at 0.6 THz. (a) Surface resis- 
tance Rs, and (b) surface reactance X$. The dashed curves are for an 
assumed 50 nm thick YBCO film. 
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Lq. The geometric inductance represents the conventional in- 
ductance of the SRR loop and can be estimated OOl to be 
Lq = 4 x 10~ n H. In contrast, the additional inductance Ls 
originates dominantly from the kinetic energy in supercon- 
ducting carriers in the YBCO SRRs. This additional SRR in- 
ductance can be calculated using the above derived YBCO 
film surface reactance X§ and by considering the geometry 
and dimensions of the YBCO SRR: L s = [(A - g) /w] {X s /to). 
Therefore, the total SRR inductance becomes L = Lq + L§. 
In order to obtain the metamaterial resonance frequency using 
Eq. we estimate the SRR capacitance C = 1.5 x 10~ 15 F, 
from the above estimated geometric inductance Lq and the 
simulated resonance frequency coo = 271 x 0.62 THz assuming 
perfect conducting SRRs (i.e. Ls = and R = ). The cal- 
culated temperature-dependent metamaterial resonance fre- 
quency is plotted in Fig.[TJc) along with the experimental and 
simulation results. Again, the theoretical result reproduces the 
frequency tuning features, though the overall tuning range is 
about half of the experimental and simulation data. 

The above calculations show that the temperature- 
dependent SRR resistance and additional inductance, due 
to the temperature-dependent complex conductivity of the 
YBCO film, play an important role in the resonance switching 
and frequency tuning. Eqs. Q and |6]) further reveal that, for 
a fixed value of the real conductivity, which is approximately 
the case in our situation, the YBCO film surface reactance, 
and therefore the additional SRR inductance, reach the maxi- 
mum value when the imaginary conductivity is approximately 
equal to the real conductivity, and vice versa. This is consis- 
tent with the experimental observations, where the metamate- 
rial resonance frequency shifts to the lowest value when the 
real and imaginary parts of the YBCO complex conductivity 
cross each other. 

The results in Fig.[2]suggest that metamaterials made from 
thinner YBCO superconducting films will have a lower res- 
onance frequency, and will be more efficient in resonance 
switching and frequency tuning. In order to verify this pre- 
diction, we fabricated and characterized a second metamate- 
rial sample from 50 nm thick YBCO film. The temperature- 
dependent transmission spectra are shown in Fig. [3] The reso- 
nance frequency at 20 K is measured to be 0.48 THz, which is 
significantly lower than that in the metamaterial sample from 
180 nm thick YBCO film. When temperature increases, the 
resonance frequency continuously shifts to lower frequencies. 
It becomes 0.31 THz at 78 K, achieving a tuning range of 
35%. We did not observe the back shifting of resonance fre- 
quency due to the high resistance at temperatures above 80 K, 
which already completely damps the metamaterial resonance. 

In conclusion, we have fabricated and characterized electric 
SRR-based metamaterials from high temperature supercon- 
ducting YBCO films. We observed temperature induced meta- 
material resonance switching and frequency tuning, which 
can be reproduced by finite-element numerical simulations us- 
ing the experimentally measured complex conductivity of the 
YBCO film. We found that both the temperature-dependent 
real and imaginary parts of the complex conductivity of the 
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FIG. 3: (color online). Temperature-dependent THz Transmission 
amplitude spectra of the 50 nm thick YBCO metamaterial. 

superconducting film have to be consistently considered in or- 
der to achieve a correct interpretation. A theoretical model 
has been developed, taking into account the SRR resistance 
and additional inductance. Our modeling calculations were in 
good agreement with experimental observations and numeri- 
cal simulations, and further predicted more efficient resonance 
switching and frequency tuning with thinner YBCO metama- 
terials, which was also verified in experiments. We expect 
that such resonance tuning in superconducting metamaterials 
could also be realized dynamically through application of op- 
tical excitation, electrical currents, and/or magnetic fields. Al- 
though high temperature superconducting metamaterials may 
not be able to essentially address the loss issue at THz fre- 
quencies and beyond, they should enable the development of 
novel, multi-functional metamaterials. 
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